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lO ■ Abstract 

O ' 

£S) ' Using the Nahm transform we investigate doubly periodic charge one SU (2) in- 

stantons with radial symmetry. Two special points where the Nahm zero modes have 
softer singularities are identified as constituent locations. To support this picture, 
the action density is computed analytically and numerically within a two dimen- 
sional slice containing the two constituents. For particular values of the parameters 
the torus can be cut in half yielding two copies of a twisted charge ^ instanton. Such 
objects comprise a single constituent. 

•i-H 

X 

Topologically non-trivial objects play a pivotal role in most confinement scenarios. A 
prominent example being the dual Meissner effect via the condensation of magnetic mono- 
poles. To date, there is no gauge independent way of identifying the monopole content of 
a given gauge potential. However, some progress has been made with regard to instantons. 
In particular, charge one SU(N) calorons, or periodic instantons, can be viewed as bound 
states of N monopole constituents [I], f2|. This identification does not hinge on a particular 
gauge choice as the constituents are clearly visible as peaks in the action density. Presum- 
ably, similar results hold for 'higher' tori, that is, instantons on T 2 x M 2 (doubly periodic 
instantons), T 3 xl and the four torus T 4 . In the latter case twists must play a role since 
untwisted charge one instantons do not exist Q. 

There is an almost complete lack of explicit results for multiply periodic instantons in 
the literature. However, the existence of higher charge instantons on T was established by 
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Taubes. More recently, doubly periodic charge one instantons have been discussed 0, |j. 
In II doubly periodic charge one SU(2) instantons with radial symmetry in the non- 
compact M 2 directions were considered. Under the Nahm transformation these instantons 

~ 2 

are mapped to abelian potentials on the dual torus T . Our basic approach is to start with 
these rather simple abelian potentials and then Nahm transform to recover the original 

SU (2) instanton. This involves solving certain Weyl-Dirac equations; for each x in T 2 x M 2 

~ 2 

one has a different Weyl equation on T . The Weyl zero modes were determined explicitly 
for a two dimensional subspace || of T 2 x M 2 (this subspace corresponds to the origin of 
M 2 ), see also §. 

Although this falls short of a complete solution, this subspace is, by virtue of the 
radial symmetry, exactly where any constituents are expected to lie. In this letter we use 
the explicit zero modes to develop a constituent picture of the associated doubly periodic 
instantons. There are two points in the subspace where the zero modes have different 
singularity profiles. These are obvious candidates for constituent locations. The behaviour 
of the field strengths in the vicinity of the proposed constituent locations is investigated 
analytically and numerically. It is possible to arrange that the two constituents are identical 
lumps. If we further fix their separation to be a half period, the torus can be cut in two. 
In each half of the torus we have a twisted instanton of topological charge | comprising a 
single constituent. 

Before we specialise to T 2 x M 2 let us briefly recall how the Nahm transform is formulated 
on T 4 [[J. Start with a self-dual anti-hermitian SU(N) potential, A^(x), on a euclidean 
four-torus with topological charge k. A gauge field on T 4 is understood to be an M 4 potential 
which is periodic (modulo gauge transformations) with respect to x^ — * x^ + L M (fi = 
0, 1, 2, 3), the L M being the four periods of the torus. The next, apparently trivial, step is 
to turn the SU(N) instanton into a U(N) instanton by adding a constant U(l) potential, 
Afj,(x) — > An(x) — iz^, where the are real numbers. We can regard the z^ as coordinates 

~ 4 

of the dual torus , T , since the shifts z^ — ► z^ -\- can be effected via periodic U{1) 

gauge transformations. 

Now consider the U(N) Weyl operator 

D Z (A) = a^(A), D*{A) = d» + A»(x) - iz», (1) 

= (l,iTi,iT 2 ,ir 3 ) where the Tj are Pauli matrices. Provided certain mathematical tech- 
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nicalities are met D\(A) = —a^D^A) has k square-integrable zero modes ifj l (x;z) with 
i — 1,2, k. The Nahm potential is defined as 



A] 3 . 



^(z)= J t /x^\x;z)^(x;z). (2) 

Here the zero modes are taken to be orthonormal. Remarkably, A(z) is a U(k) instanton 
on the dual torus with topological charge N. If we execute a second Nahm transformation 
on A, the original SU (N) instanton will be recovered. 

Formally, one can obtain the T 2 x M 2 Nahm transform by taking two of the periods, say 
L and L 3 , to be infinite. Given an SU (N) instanton periodic with respect to x 1 — > x 1 + L 1 , 
x 2 — * X2 + L 2 its Nahm transform is 

A»(z) = f d 4 x^(x;z)^(x;z) p = 1, 2 

A$(z) = / d 4 x^(x;z)ix tM i) j (x;z) fj, = 0, 3 (3) 

where the ip l (x; z) with % = 1, k are orthonormal zero modes of D\(A). Note that we 
may gauge z and z 3 to zero. The Nahm potential is a self-dual U (k) gauge field with N 

~ 2 

singularities in T (periods 2ir / 'L\ and 2tt/L 2 ). 

One can regard Ai(z) and A 2 (z) as the components of a two dimensional gauge poten- 
tial, and combine A (z) and A 3 (z) into a 'Higgs' field $(z) = \(Aq(z) — iA 3 (z)). The next 
step is to seek solutions of the dimensionally reduced self-duality (or Hitchin) equations. 
In the one-instanton sector this is rather straightforward since the corresponding Nahm 
potential is abelian. What is more tricky, however, is to execute the second Nahm trans- 
formation to recover the instanton itself. We shall restrict ourselves to the special case of a 
zero Higgs field. This means that the corresponding instanton will be radially symmetric: 
local gauge invariants such as the action density depend on x\, x 2 and r = a/xq + x\ only. 
With this restriction the self-duality equations are just 

Fyy = OyAy ~ G )y Ay = 0. (4) 

Here we have used complex coordinates; y = Z\ + iz 2 , y = Z\ — iz 2 , d y = \{d zi — id Z2 ), 

dy = \{d Zl + id Z2 ), A y = \(A 1 (z) - iA 2 (z)) and Ay = \{A 1 (z) + iA 2 (z)). Consider the 
ansatz 

Ay = dy(f), Ay = ~dy(/), (5) 



which gives F yy - = —2d y dy(f). Then, in order to satisfy (f|), (f> must be harmonic except at 
two singularities, since we are aiming for an SU(2) instanton. A suitable <fi satisfies 



K + d 2 Z2 ) fa) = -2nn [6 2 (z - u>) - 5\z + uj)} , 



(6) 



where k is a constant and ±u are the positions of the two singularities (we have used 
translational invariance to shift the 'centre of gravity' of the singularities to the origin). 
The delta functions should be read as periodic (with respect to Z\ — > Z\ + 2n/L 1 and 
z 2 z 2 + 27[/L 2 ). Physically, the Nahm potential describes two Aharonov-Bohm fluxes 
of strength k and — k threading the dual torus. They must have equal and opposite 
strength to ensure a periodic A. We may assume that k lies between and 1 since it is 
possible via a (singular) gauge transformation to shift k by an integer amount (under such 
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a transformation the total flux through T remains zero). 

One can write explicitly in terms of Jacobi theta functions 
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The theta function is defined as 
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and has the periodicity properties 9(w + l, r) = 9(w, r) and 9(w + r, r) = e 
In each cell 9(w,r) has a single zero located at the centre of the torus (w 
have chosen the constant term in (0) so that the integral of (f) over the dual torus is zero. 
This renders <f>(z) an odd function, (j){—z) = —<p(z). The Nahm potential derived from ([!]) 
was obtained in H via the ADHM construction. Here the flux strength k is related to the 



ADHM 'size', A, of an instanton centred at x M = 0, 
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Although we do not directly use the ADHM formalism in the present paper the relation 
(|) proves useful in interpreting our results. 

The SU(2) instanton we seek, A^x), is the Nahm transform of the abelian potential 

A(z) 



/ 
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d 2 z^(z;x)—i/j q (z;x) 
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where the ip p (z; x),p — 1, 2 are orthonormal zero modes of 
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In addition to the complex coordinates y, y on T we have introduced two sets of complex 
coordinates for T x M 2 ; in the 'parallel' directions x\\ = x% + 2x2, X|| = xi — 2x2, and in 
the 'transverse' non-compact directions x± = xo + 223, = xo — 2x3. 

When x_i_ = the Weyl equation decouples and the two zero modes have a simple form 
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ip\z;x) 







ip 2 (z;x) 




e-+toG + (z-w) 
where G±(z) are periodic Green's functions satisfying 

(-id y - G + (z) = \5\z), {-idy - \x\\) GL(z) = \8 2 {z). 

G-(z) has a theta function representation 
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with 6'{w,t) = d w 6(w,r). The corresponding result for G + (z) can be obtained via 
G + (z) = G*_(—z). The zero modes have square-integrable singularities at both fluxes: 
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Table 1: Singularities of the x± = zero modes 

There are two values of xm (in each copy of T 2 ) for which Tab. 1 does not hold. When 
x\\ = the Green's function's G± do not exist. This does not mean there are no zero 
modes, indeed one can see that 



1)\z;x = 0) 
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ip 2 (z;x = 0) 
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are solutions of the x\\ = Weyl equation, ip 1 has the expected square-integrable singularity 
at z = — u, but for z = u, ip (z;x = 0) is zero. On the other hand ip (z,x = 0) diverges 
at z = uj but not at z = —uj. The finiteness of G-{2uj) was used to derive Tab. 1. But if 
it is zero, ip l will lose its singularity at z = —uj. From the theta function representation 
of G-(z) one can see that there is exactly one value of xu in T 2 where this occurs, wxu = 
—iL\L 2 {uj\ + 1002). We wish to investigate whether these 'soft' points are constituent 
locations. 

Notice that the singularity profiles of ip l and ip 2 are exchanged under the replacement 
k — > 1 — k. This suggests that the constituents are exchanged under this mapping. That 
is, if there are indeed lumps at the two points, then k —>■ 1 — k swaps the two lumps. The 
following result formalises this idea 

F llv {x\\,x±, k) = V~ 1 (x)F fW (~x\\ + ^r^-(ui + iuj 2 ), -x±, 1 - kJ V(x), (16) 

where V(x) is some U(2) gauge transformation. The proof of Eq. (O) goes as follows: 
make the change of variables z — > — z in (TO). The zero mode ip p (—z; x) satisfies the same 
Weyl equation as ip p (z; x) except that the signs of k and the x M are flipped. Under periodic 
gauge transformations A(z, —k) and A(z, 1 — k) are equivalent up to a constant potential 

A^z, -«) = U-\z) (d, + A^z, 1 - «) + B^j U(z), (17) 

where ~kB\ = —iLiL 2 u)2 and nB 2 = iL\L 2 LO\. In the Weyl equation B^ can be absorbed 
into x\ and x 2 . Thus U^fyipv^z, —x\\ — iLiL 2 {uji + iu 2 )/n, —x±) satisfies the same Weyl 
equation as xjj p (— z; x\\, x±). 

We have established that any lumps centred at the constituent points have a simple 
exchange property. To actually see the lumps we must compute the field strengths. We 
can do this explicitly in the x± = slice, since here we have the exact zero modes (|T2|). 
Before we insert these modes into (|10D, they must be normalised. This can be done by 
dividing both modes by ^/p, where 

p = [ dhe^lG^z + co)] 2 . (18) 

In fact, one can write two components of the gauge potential solely in terms of this nor- 
malisation factor 

4b N = -| r 3«9x|| log/?, As u = §r 3 <9 S|| logp, (19) 
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so that 

^11*11 = r 3 d xn dx u hgp. (20) 

The computation of the other components is more involved since the derivatives of the zero 
modes with respect to x± are required. However, the field strengths can be more easily 
accessed via Green's function techniques. Here we just quote the results; more details will 
be given elsewhere 

F^ s± = 2m(n + iT 2 )Kp^. (21) 

where 

v= f d 2 zG + (iu - z)e 2 ^ z) G„(z + u). (22) 
Jf 2 

The other components are fixed by self-duality, i.e. F x±x± +F X .. X .. = and F X .. XJ _ = 0. Note 
that p is dimensionless, real and periodic (with respect to x\ — > X\ + L%, x 2 — > x 2 + L 2 ), 
while v is dimensionless, complex and periodic up to constant phases. Both p and v 
diverge at Xu = 0, but the field strengths should not. A careful analysis of p and v in 
the neighbourhood of X|| = shows that the field strengths are well defined at this point. 
Alternatively, one can just note that p and v are well behaved at the second soft point and 
invoke (Pp. 

Since we have all field strengths within the x± = slice we can compute the action 
density 

x 2 . „ 2 

(23) 



- ±Tri^F^ = 16 (d Xu d X]} \og P y + 1287rVp 2 
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We have provided integral representations of p and v where the integrands are expressible in 
terms of standard functions. Using these results, we have made numerical plots of the action 
density for various values of k and u, and for the periods we have taken L\ = L 2 = 2ir, i.e. 
a square torus. The most symmetric and ui = uj 2 = \ is plotted in Fig. |l]. 

Two equal sized lumps are observed, one at xu = 0, the other at xu = (1 + Decreasing 
k at fixed uj reduces the size of the x\\ = constituent but increases its contribution to the 
action. Conversely, the second constituent becomes larger but contributes less to the action 
density. These two effects tend to (rather quickly) flatten the second peak as k is reduced. 
Here we provide plots with k = ^ (Fig. |^) and k — | (Fig. |), again with u\ — uj 2 — |. 
Already at k = | we see that the x\\ = (1 + i)ir constituent is more of a plateau than a 
peak. If k is decreased much further the first peak will dominate completely; this peak 
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will be essentially that of a BPST instanton on R 4 with a scale parameter given by the 
ADHM formula @. Indeed, if we take (H) at face value, and use (|16|), our two constituents 
are BPST instanton cores, one with scale A = a/ k,LiL 2 /tt at x\\ = x± = 0, the other with 
scale A = a/(1 — l^LiL^pK at the second soft point. When L\ = L 2 at least one A is of 
the order of the size of the torus in line with the strong overlap seen in the plots. 

The peaks are best resolved at k — |, when the constituents are identical. The overlap 
can be avoided by taking one period, say L 2 , smaller than the other, and then taking a 
sufficiently large constituent separation in the x\ direction. However, then the size of one 
or both cores will exceed L 2 , tending to suppress or smooth out the dependence on x 2 . In 
this case the constituents become periodic monopoles. Actually, it is possible to see some 
hint of this monopole limit even for a square torus; taking u\ = |, lo 2 = 0, k = ~ with 
L\ = L 2 = 2ir the two constituents (at xu = and xu = in) overlap so strongly that they 
merge into a monopole worldline, in that the dependence of the field strengths on x 2 is 
very weak. 

The has another interesting feature. If we choose the constituent locations 

so that they are separated by half periods the charge one instanton can be 'cut' to yield 
a charge | instanton (see also 0). This happens when (lji,uj 2 ) is (^tt/Lx,0), (0,^ti/L 2 ) 
or (~7r/Li, \ti/L 2 ); the former and the latter are shown in Fig. |J After cutting we have 
a twist Z\2 = — 1 in the half torus. But to get a half-integer topological charge we also 
require a twist in the non-compact directions, = — 1, since non-orthogonal twists are 
a prerequisite for fractional instanton number. Far away from x± = 0, the potential must 
be a pure gauge 

A^x) — V- 1 {x 1 ,x 2 ,6)d fl V(x ll x 2l 6) r^oo, (24) 

where x± = re ld . The non-compact twist, Z03 = — 1, translates into a double- valued gauge 
function, V(xi,x 2 ,9 + 2n) = —V(xi,x 2 ,9). However, the explicit results we have given 
are restricted to the r = slice. What can we say about the full r ^ potential? It can 
be argued that as for r = (as well as £77(2) calorons) two functions, one real the other 
complex, are sufficient 

r 3 v* 

A x\\ = "TT^H logP - 27l7(ri - iT 2 )Kpds: ± —, 

A x± = -^-d x± \ogp + 2ni(T 1 -iT 2 )Kpd Sl —. (25) 
A p 
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Here p and v depend on r. The p and v respectively given in (|i~8f) and (^) should be 
understood as the zeroth order terms of power series in r 2 for the 'full' p and v. Assuming 
that v falls off exponentially for large r, a decay p ~ r _1 , would account for the twist. 

Doubly periodic charge \ instantons have also been found in simulations p, [|. Many 
qualitative features reported are in accord with our findings. In particular, we agree that 
these objects are single lumps with radial symmetry, and that the action density is never 
zero in the x± = plane. Another striking observation in ||, |9| is that the action density 
decays exponentially for large r. To make a comparison here, more information about the 
large r properties of p and v is required. 
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Figure 2: Plot of action density for x± = and k = X and u% = u>2 
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Figure 3: Plot of action density for x±_ = and k = §, u>i = o> 2 = \- 
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Figure 4: When k = |, the choices cji = |7r/Li, co> 2 = and tui = ^ir/Li, u 2 = |vr/L 2 
yield constituent separations allowing the torus to be cut (dashed line) to yield charge | 
instantons. 
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